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Abstract: 

An investigation of the hadronic final state in diffractive and non-diffractive deep-inelastic 
electron-proton scattering at HERA is presented, where difl^ractive data are selected experi- 
mentally by demanding a large gap in pseudo-rapidity around the proton remnant direction. 
The transverse energy flow in the hadronic final state is evaluated using a set of estimators 
which quantify topological properties. Using available Monte Carlo QCD calculations, it is 
demonstrated that the final state in diffractive DIS exhibits the features expected if the in- 
teraction is interpreted as the scattering of an electron off a current quark with associated 
effects of perturbative QCD. A model in which deep-inelastic diffraction is taken to be the 
exchange of a pomeron with partonic structure is found to reproduce the measurements well. 
Models for deep-inelastic ep scattering, in which a sizeable diffractive contribution is present 
because of non-perturbative effects in the production of the hadronic final state, reproduce 
the general tendencies of the data but in all give a worse description. 
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1 Introduction 



In high energy physics, the parton model together with the theory of strong interaction, Quantum- 
Chromo-Dynamics (QCD), have been shown to provide a good description of a variety of different 
processes involving hadrons in the final and/or initial state (e.g. in e^e~' collisions Q, in hadron-hadron 
collisions!^ ^^'^ lepton-hadron scattering [Q). However, so far QCD has not been able to make predic- 
tions, derived from first principles, for a class of hadron-hadron collisions known as diffractive scattering. 
This class is characterised experimentally as peripheral collisions between the incoming hadrons, whereby 
the scattered hadrons keep their original identity or break up dominantly into a system of low invariant 
mass. This generally leads to experimentally observable large rapidity gaps in these collisions. Regge 
theory applied to hadronic interactions models diffractive collisions [Q as the exchange of an object called 
the pomeron (IP), which carries only the quantum numbers of the vacuum and thus no colour charge. 
Attempts have been made to study the nature of diffractive exchange, and it has been suggested that 
the pomeron has partonic structure indications for which were found in proton-antiproton collisions 
by the U A8 experiment Q . 

At the turn-on of the electron-proton (ep) collider HERA, at DESY, Hamburg, a class of deep- 
inelastic scattering (DIS) events was observed which exhibited in the hadronic final state an unusu- 
ally large rapidity gap (LRG) with almost no hadronic energy flow around the direction of the proton 
remnant]^, |j. Recent cross-section measurements |l^ have shown that these events exhibit charac- 
teristics similar to those of diffractive hadronic collisions. At low Bjorken-a; (see below for kinematic 
definitions) the diffractive contribution to the total DIS cross-section, and thus to the proton structure 
function F2, has been quantified in the form of a diffractive structure function Fj^^"^"*^, ^]. The de- 
pendence of i^2^^^-' on the appropriate deep-inelastic kinematic variables xp, Q^, and P (see below) has 
been measured. It demonstrates both that the diffractive deep-inelastic ep process can be interpreted 
as deep-inelastic scattering of the incident electron off a colourless object coupling to the proton in the 
diffractive ep interaction, and that the structure of this object is consistent with being of a partonic 
nature. Thus deep-inelastic scattering may be pictured as in figure □ (a), and the diffractive contri- 




Flgure 1: Schematic picture for deep— inelastic scattering: non-diffractive (a) and diffractive (b). 

bution to it as in figure |l] (b), where the process is considered in the phenomenology which is used to 
quantify diffractive interactions in high energy hadron-hadron physics. In each case a virtual photon (in 
the kinematic region investigated, contributions from exchange can be neglected) probes a hadronic 
object. For non-diffractive DIS this object is a proton, whereas in the diffractive case this is taken to be 
a colour neutral object emitted from the proton. The proton remnant system in the latter case remains 
colourless and thus a gap in rapidity is produced. Note that experimentally we observe events which 
have a rapidity gap, as defined below. Throughout this paper we will however use the terms "diffractive" 
event and "LRG event" synonymously. 



Using the four-momenta k of the incident electron, q of the virtual photon and P of the incident 
proton the following kincmatical variables can be defined: 

Q' = -l', ^=S-q^ W'^iP + qf, (1) 

where is the squared virtuality of the photon and x is the Bjorken scaling variable, which in the naive 
Quark-Parton Model (QPM) can be interpreted as the fraction of the proton's momentum carried by the 
struck quark. In the rest frame of the proton y is the fraction of the electron energy transferred to the 
proton. W is the invariant mass of the photon-proton system and is equal to the invariant mass of the 
total final state excluding the scattered electron. 

Additional kinematical variables can be defined using the four-momentum P' of the colourless remnant 
(either a nucleon or a higher mass baryon excitation) in the final state: 

xjp p= ^ and t={P-P'f. (2) 

q-P ^ 2q-{P-P') ^ ' ^ ' 

The variables x, xp and (3 are related via x — xpP. Defining Mx to be the invariant mass of the virtual 
photon-pomeron system (i.e the invariant mass of that part of the final state not associated with the 
colourless remnant and separated from the latter by a LRG), xp and (3 can be written as 

+ Ml - 1 

where Mp is the mass of the proton. When and \t\ are small (M^ < and < W'^ and \t\ ^ 
and < Ml), Xp can be interpreted as the fraction of the proton's four-momentum transferred to the 
pomeron, and (3 can be viewed as the fraction of the pomeron's four-momentum carried by the quark 
entering the hard scattering. In the kinematic region under investigation both and \t\ can be neglected 
and therefore xp and (3 can be calculated from Af^, and as 

The consistency of the dependence of F^*-"^^ on (3 and with a partonic interpretation of the 
pomeron [|[ |l^ implies that the hadronic final state in deep-inelastic diffractive scattering is expected 
to show evidence for parton production and effects of perturbative QCD of a nature similar to deep- 
inelastic ep scattering at appropriate {x ~ (3) Bjorken-a;. This picture can be further tested through 
detailed measurements of the hadronic final state. In this paper we present an analysis of the hadronic 
final states for diffractive and non-diffractive DIS, and compare the results with QCD inspired models. 
The analysis concentrates on the observed energy flow in the laboratory frame relative to the expected 
direction of the struck quark. In the naive Quark-Parton Model this direction can be calculated from 
the measurement of the scattered electron using four-momentum conservation. Such prediction of the 
quark direction is not possible in e+e^and hadron-hadron collisions. There the measured hadronic flnal 
state must itself be used, e.g. by defining jet directions which then can be related to parton directions. 
In|pT[ the ZEUS collaboration compared the final state of events with a LRG to that of events without a 
gap using energy flow measurements and concluded that in DIS with a LRG QCD radiation is strongly 
suppressed. Recently a comparison of charged particle spectra for the two classes of events was performed 
by the ZEUS coUaborationp^ , where they observed similarities between DIS with a LRG at HERA and 
DIS at lower W (w Mx). 

In this analysis, properties of the observed energy flow are defined which are sensitive to effects of 
perturbative QCD. The evolution of these properties with the kincmatical variables is investigated. The 
data presented are corrected for detector effects. From the diagrams shown in figure ^ it is expected that 



the different fcinematics in the two cases impiy differences in the avaiiable phase space for QCD effects 
in the hadronic final state. The invariant mass of the system built from the photon and the probed 
object is given by W for non-diffractive DIS and Mx in the difTractive case, where Mx «C W. The 
relevance of this scale for hadron production in the diffractive process will be tested by considering DIS 
at a reduced proton beam energy. This feature is also a key ingredient of several more sophisticated 
models for diffractive DIS, such as RAPGAP[p3[. Here it is assumed that the hard scattering occurs off 
a partonic constituent of a pomeron emitted from the proton as shown in diagram (b). QCD effects 
of parton showers and hadronization are included in this model. However, alternative explanations for 
the diffractive process exist. The two models LEPTO[^ and HERWIG[p^, originally developed for 
non-diffractive DIS, in the most recent versions also generate events with a leading colourless remnant. 
These models do not involve explicitly the emission of a pomeron but instead produce the gap through 
non-perturbative QCD effects in the evolution of the final state produced by deep-inelastic scattering of 
a partonic constituent of the proton. The predictions of these models will be compared with the data. 

The paper is organized as follows. After a short description of the HI detector (section 2), the data 
taking and the event selection are briefly discussed (section 3). Then are described the models of the 
hadronic final state used in deep-inelastic scattering (section 4), the procedure to correct for detector 
effects, and the sources of systematic errors (section 5). In section 6 the characteristic properties of the 
energy flow are explained and the estimators used are defined. Results and the comparison with different 
model predictions are presented in section 7. The paper is summarized in section 8. In an appendix the 
measured values of the estimators are listed together with their statistical and systematic errors. 

2 The HI Detector 

The general layout of the HI detector (described in more detail in is as follows: the interaction region 
is surrounded with tracking devices which measure the momenta of charged particles and reconstruct the 
position of the interaction point. These tracking detectors are enclosed by a calorimetric region which 
allows the measurement of the energy and the direction of charged and neutral particles. All these 
detectors are situated within a magnetic field of 1.15 T, generated by a superconducting coil. The flux 
return yoke of this coil is instrumented to identify muons and to measure energy escaping the main 
calorimeters. In the following the components of the HI detector of particular relevance to this analysis 
are briefly presented. In the coordinate system used, the proton beam direction defines the +z-axis. 
The region of polar angle O<0<7r/2is called "forward region" and corresponds to positive values of 
pseudo-rapidity — ~ In tan |, whereas ?/ < in the "backward region" {n/2 < 9 < n). 

The measurement of charged particle tracks and the determination of the interaction point is made 
with a system of interleaved drift and multiwire proportional chambers covering the central and forward 
regions of the detector (7° < 9 < 165° corresponding to —2.0 < rj < 2.8) over the full azimuthal 
range. The backward region is equipped with a proportional chamber (BPC) measuring charged particle 
positions in the angular range 155° < 9 < 174.5° in front of the backward calorimeter (BEMC). 

The main calorimeter is a fine-grained liquid argon (LAr) calorimeter , covering a range in polar 
angle from 4° to 155° (—1.51 < rj < 3.35). It consists of an electromagnetic section with lead absorbers 
(20 - 30 Xq) and a hadronic part with steel absorbers, giving a total interaction length of 5 - 8 A. The 
energy resolution as measured in test beams JTsI is cr/E w 12%/y/E for electrons and w for 
hadrons, with E in GeV. The energy scale is known to about 2% for electrons and to about 5% for 
hadrons. 

Covering the 77 range —3.35 to —1.35 (151° < 9 < 176°), the backward electromagnetic calorime- 
ter (BEMC) allows the measurement of the scattered electron for low DIS events (5 GeV^ < < 



100 GcV^) and provides information on hadrons scattered in this region. The BEMC is a lead-scintillator 
sandwich calorimeter (22.5 Xq or 1 hadronic absorption length), with a resolution of s:! 10%/^/E for elec- 
trons and an uncertainty in the energy scale for electrons of about 1.7%. Located behind the calorimeter 
is a time-of-flight system, consisting of scintillators which veto upstream interactions of the proton beam. 

The very forward region is equipped with a copper-silicon calorimeter (PLUG), covering a range in rj 
from 3.54 to 5.08. The longitudinal thickness is 4.25 hadronic absorption lengths. A muon spectrometer 
surrounds the beam pipe outside the flux return yoke in the forward direction. It consists of a toroidal 
magnet sandwiched between two sets of drift chambers and is used to detect charged particles from 
ep interactions in the range 5.0 < i] < 6.6 by means of secondary particles which they produce from 
interactions in the beam pipe and adjunct material. The latter two detectors are used in this analysis to 
tag particle production in the forward region close to the proton beam direction. 

For the determination of the luminosity the process ep — » ep7 is used, where the electrons and photons 
are measured in two calorimeters located far downstream of the detector in the electron beam direction 
(the electron "tagger" at ^ = —33 m and the photon "tagger" at z = —102 m). 

3 Data Taking and Event Selection 

In 1993 the HERA collider operated with a 26.7 GeV electron beam and an 820 GeV proton beam, giving 
a centre of mass energy of y/s = 296 GeV. The data used for this analysis correspond to an integrated 
luminosity of 294 nb^^. The events used here were triggered by requiring a cluster with an energy larger 
than 4 GeV in the BEMC and no veto from the time-of-flight system. This trigger has an efficiency of 
99% for events containing an electron with an energy of more than 10 GeV in the angular acceptance of 
the BEMC and provides a sample of DIS events at low (5 QeV^ < < 100 GeV^). 

The measurement of the energy E'^ and the polar angle 9e of the scattered electron is used to determine 
the kinematical variables = iEeE'^ cos^(0e/2) and y — 1 — (E'^/Ee) sin^(0e/2), where Ef, is the electron 
beam energy. Bjorken x and the square of the invariant mass of the total hadronic system, W'^, can be 
calculated from and y using the relations x = Q^/{ys) and = sy — Q^. 

The selection of deep-inelastic scattering events follows closely that used in the recent measurement of 
the deep-inelastic structure function ^2(1, (5^)|l^ and the diffractive structure function F^''^\x, Q'^^xp) 
by HI. The main requirements are an electron candidate in the backward calorimeter, a reconstructed 
interaction point, and a minimal invariant mass of the hadronic final state: 

• The electron candidate is required to have an energy E'^ > 10.6 GeV and to be found within the 
angular acceptance of the BEMC: 155° < 6e < 173°. In addition, electron identification cuts to 
suppress background from photoproduction events arc made using the information obtained from 
the cluster radius and the matching of the cluster with a charged particle signal in the BPC. 

• A reconstructed vertex close to the nominal interaction point is demanded: \zvertex — ^"e™e""'l < 
30 cm. 

• To ensure accurate reconstruction of the kinematics from the detected scattered electron, y > 0.05 
is required since the reconstruction deteriorates at low values of y. 

The selected events cover the kinematical range 10~^ < x < 10^^ and 7.5 GeV^ < < 100 GeV^ 
with an average value for of 23000 GeV^ (4300 GeV^ <W^ < 53000 GeV^). 



For the measurement of energy flow and the invariant mass of the final state, clusters reconstructed 
in the LAr calorimeter and in the BEMC were used. In the BEMC, only clusters with an energy larger 
than 400 MeV are considered. 

The sample of DIS events obtained contains events with and without a large gap in rapidity. DIS 
events without a gap are selected by demanding 

• a minimal energy deposit in the forward region measured in the LAr calorimeter: Eforward > 
0.5 GeV (as used in |2l)), where Eforward is the summed energy in the the region 4.4° < 9 < 15°, 
corresponding to 2.03 < i] < 3.26. 

This requirement together with the DIS selection described above results in 15242 events. The only signi- 
ficant background to these events from other ep interactions is due to photoproduction which contributes 
about 9% of the events at the lowest value of x « 2 • 10^^ and can be ignored for values of x > 4 • 10^^. 

For the selection of DIS events with a large rapidity gap, as described in detail in |9| , the existence 
of a region around the proton remnant direction with almost no hadronic energy flow is required. The 
detector components used for this selection give access to the very forward region (up to a pseudo-rapidity 
of ry w 6.6): 

• The energy deposited in the Plug calorimeter has to be smaller than 1 GeV and the number of 
reconstructed hit pairs in the forward muon spectrometer has to be smaller than 2. In addition 
Vmax,LAr < 3.2 is required, where rjmax,LAr is the maximum pseudo-rapidity of all clusters in the 
LAr calorimeter with E > 0.4 GeV. 

Applying these cuts together with the DIS kinematical cuts leads to a sample of 1721 events with a LRG. 
This sample consists of events where a leading hadron or cluster of hadrons Rc in the forward direction 
escapes detection by remaining in the beam pipe. The acceptance in the invariant mass Mj^ of Rc is 
specified by the forward detectors used to define the sample. As Mn increases from the mass of the 
proton to 4 GeV the acceptance decreases from 100% to less than 5%; above 4 GeV the acceptance is 
always less than 5%. In the central detector the remainder of the final state, separated by a gap from 
Rc, is detected. The invariant mass Mx of this system defines through equation 4 the kinematic variable 
xp for each event. The data sample^ used covers the range 2 • 10^* < xp < 2 • 10^^. In all results 
presented, diffractive DIS is taken to be for x < xp < 0.02. 

4 Monte Carlo Models for the Hadronic Final State 

The Monte Carlo models for DIS can be separated into three parts. The hard interaction of a virtual 
boson with a parton is simulated using the leading order electroweak cross-section for parton scattering 
including the first order QCD correction given by the exact 0(as) matrix elements (these are: the Born 
term for boson-quark scattering, hard gluon Bremsstrahlung, and the boson-gluon fusion process). A 
parton cascade includes higher order QCD corrections to generate additional partons. The resulting 
coloured partons hadronize to give the observable particles. 

The models differ mainly in the details of the parton cascade and the phenomenological description 
of the hadronization. The following description of the models emphasizes these different approaches. 

^The sample of events with a gap and that without a gap are not completely disjunct. About 5% of events without a 
gap are also classified as events with a gap. This however does not effect the conclusions drawn. 



LEPTO 14 1 uses the the leading-log approximation based on the Altarelli-Parisi evolution equa- 
tions 1^ for the parton showers. The fragmentation is done via the Lund string model||2^ as 
implemented in JETSETQ. 

The new version (6.3) of LEPTO differs from the previous version (6.1) mainly in two aspects. 
Firstly, the treatment of scattering involving a sea-quark has been modified, motivated by the poor 
description of the measured transverse energy flow in the forward region|^l| given by version 6.1. 
Secondly, the possibility of colour rearrangement in the final state through soft gluon exchange 
with negligible change in the momenta of the partons has been introduced, allowing the generation 
of events with a leading colourless remnant. This soft colour interaction is a non-perturbative 
interaction of the coloured quarks and gluons with the colour medium of the proton. This way 
of generating a colour-neutral subsystem is similar to the ideas of Buchmueller and Hebecker|p^, 
who performed a calculation of the diffractive cross-section on the parton level. The size of the 
diffractive contribution is determined by a probability^ that such a colour exchange occurs between 
two colour charges, leading for part of the cross-section to the formation of colour singlet subsystems 
separated in rapidity. 

ARIADNE is a generator for QCD cascades only. In this analysis, two versions (4.03 and 4.07) 
are used. The former version is used to calculate the correction for detector effects for events 
without a gap (see next section). Previous analyses j2^, ^ have shown that this version gives a 
reasonable description of the data. The version 4.07 is used in comparison with non-diffractive 
DIS. For the modelling of the electroweak interaction the corresponding parts of LEPTO[0 are 
used. Gluon radiation is performed in ARIADNE by an implementation of the Colour Dipole 
Model [p9| . In this model, gluon emission from a quark-antiquark pair is treated as radiation from 
a colour dipole formed by this pair. In contrast to the quarks in e+e~-annihilation, the proton 
remnant in DIS is assumed to be extended. This leads to a suppression of the phase space for 
gluon radiation. In version 4.07 the struck quark is considered to be extended as well. The latter 
modification was motivated by the disagreement between HERA data and the ARIADNE (version 
4.03) prediction in the region of the struck quark pl|, |2^. The hadronization is done using the Lund 
string model (JETSET). In version 4.07 diffractive DIS is modelled via the emission of a pomeron 
from the proton and subsequent hard scattering on a partonic constituent of the pomeron. This 
mechanism follows the spirit of RAPGAP (see below) and is not investigated in this paper. In[^ 
Lonnblad demonstrates that in the framework of the Colour Dipole Model colour reconnections 
cannot reproduce the diffractive contribution to DIS in contrast to the LEPTO model. 

HERWIG is a general purpose generator for high energy hadronic processes. The parton shower 
algorithm (in the leading logarithmic approximation) takes into account colour coherence as well 
as soft gluon interference. The hadronization in HERWIG is performed using the concept of cluster 
fragmentation, where gluons are split non-perturbatively into quark-antiquark pairs. The latter 
are combined into colour singlet clusters, which are split further until a minimum value of the 
cluster mass is reached. In the most recent version (5.8d)|3^, as used in this analysis, HERWIG can 
generate a diffractive contribution to DIS, the size of which is sensitive to a parameter^ determining 
the mass distribution for the cluster splitting in the hadronization. 



For all the above models the MRS(H)^^ set of parton distributions for the proton was used. These 
were determined using data from various experiments including the F2 measurements from HERA made 

^The default value for PARL(7) is 0.2, this leads to a diffractive contribution of about 7% to the deep— inelastic cross 
section in the kinematic range considered. A variation of this parameter between 0.1 and 0.5 gives a change in this fraction 
between 5% and 9% but no significant changes in the measured energy flow are observed within the two event classes. 

''The chosen value for PSPLT is 0.7, this leads to a diffractive contribution of 6%. A variation of this parameter between 
0.5 and 1.0 gives contributions of 2% and 8%. However the change in PSPLT results in significant changes in the predicted 
energy flow properties for both classes of events in contrast to LEPTO. 



in 1992|3^, |3J]. In addition to the kincmatical selection as described above a cut on the summed energy 
Eforward of all particles produced in the range 2.03 < rj < 3.26 is performed by demanding E forward > 
0.5 GeV for the distributions obtained from the models in the non-diffractive case. 

It should be stressed that the LEPTO model as well as the HERWIG model in the most recent 
versions can generate a difFractive contribution to the DIS cross-section compatible with the measurement 
without explicitly involving the concept of deep-inelastic electron-pomeron scattering. These events are 
selected by demanding that rj^J^^ < 3.2, where rj^ax the maximum pseudo-rapidity of all particles with 
E > 0.4 GeV and rj < 6.6. Both models broadly reproduce the xp dependence of F^^^"^ as measured 
in|^ and thus the one of the deep-inelastic diffractive cross-section. Most of the parameters in these 
models which influence the description of the hadronic final state are restricted by measurements of non- 
diffractivc DIS, leaving at present only a few parameters free for the modelling of diffractive DIS. Events 
with a LRG also occur in the previous versions of the models but at a very small rate (being due to 
extreme fluctuations in the fragmentation). 

RAPGAP models diffractive DIS by the emission of a pomeron from the proton, described by a 
flux factor depending only on a: ip and t. The pomeron is taken to be an object with a partonic 
structure described by parton densities, which depend on f3 and Q^. This hypothesis is consistent 
with the recent measurements of the diffractive structure function F^'^'*!^, In the version 

used (1.3) the parton content of the pomeron can be chosen to be either a quark-antiquark pair or 
two gluons. Within this analysis, the parton densities p{z) are chosen to be "hard" distributions 
{[z ■ p{z)] ^ z ■ {1 ~ z)), where z is the fraction of the momentum of the parton relative to the 
pomeron. A soft parton density ([z • p{z)] (1 — z)^) is excluded by the measurement of f^^^^ 
as described in A mixture of "hard quark" and "hard gluon" densities was used such that an 
equal number of events are generated for both densities in the kinematic range considered. This 
leads to a sample in which, for /3 < 0.1, more than 70% (and for (3 > 0.9 less than 5%) of the 
events are of the "hard gluon" type. The mixture has been chosen to get a good description of the 
measured diffractive structure function F^^^^ and was found to be able to describe the measured 
energy flow. Additional partons due to QCD radiation are generated using the Colour Dipole model, 
the subsequent hadronization is performed with JETSET. The version of RAPGAP used does not 
include the evolution of parton densities with Q^. 



5 Correction for Detector Effects 

The correction of the measured energy flow for detector effects in non-diffractive DIS is done using events 
generated with the ARIADNE 4.03 127| model, which have been passed through a simulation of the HI 
detector response. They are reconstructed in the same way as is done for the data. For diffractive DIS 
the same procedure is apphed using the RAPGAP 1.3||l^ model. 

For each distribution shown, a set of bin-by-bin correction factors is calculated by forming the ratio 
of the bin-contents for the generated events (at the particle level) to the corresponding bin-contents for 
the reconstructed events (from the detector simulation, i.e. at the detector level). To obtain the corrected 
value, the raw data bin-contents have to be multiplied by the correction factor. The derived factors vary 
only moderately from bin to bin and have values typically between 0.8 and 1.2. 

In the determination of the systematic error associated with the correction factors, the following 
effects were considered: 

• Variation of the energy scale for the scattered electron in the BEMC. The energy scale is known to 
±1.7%. 



• Uncertainty in the polar angle of the scattered electron. An error of ±2 mrad was taken into 
account. 



• Uncertainty in the hadronic energy scale in the LAr calorimeter. From studies using pT^balance 
between the scattered electron and the hadronic final state, the energy scale is known to an accuracy 



• Uncertainty in the hadronic energy scale in the BEMC. The hadronic energy scale in the BEMC 
was assumed to be known to ±20%. 

• Dependence on the model used for corrections. For events without a LRG, a comparison between 
the correction factors obtained using the ARIADNE and the LEPTO model was performed. In the 
case of DIS with a LRG, the RAPGAP model assuming either a quark parton density or a gluon 
density alone was studied. 

• The effect of initial state photon radiation off the electron has been estimated with the DJANGO 



• The effect of background from photoproduction {Q^ « 0) has been investigated using HI photopro- 
duction data in which the scattered electron is detected in the electron calorimeter of the luminosity 
system and an electron candidate is found in the BEMC. The values of the estimators obtained were 
found to be smaller than those in the DIS data. Using a value of 9% for the contamination from 
photoproduction at the lowest value of a; « 2 • 10~'*|jl^ (contamination negligible for a; > 3 • 10^"*) 
an asymmetric contribution to the systematic error is obtained. 

All contributions have been added quadratically to give a value of the systematic error for each bin 
considered. The error bars shown contain the statistical error (inner bars) as well as the total error (full 
error bar) which has been obtained by adding statistical and systematic errors in quadrature. 

6 Characteristic Properties of the Energy Flow 

Within the framework of the naive Quark-Parton Model, the measurement of the four-momentum of the 
scattered electron determines the direction of the quark struck in the deep-inelastic scattering. Using 
conservation of four-momentum and assuming that the partons are massless and the proton remnant 
has negligible transverse momentum, the polar angle of the struck quark can be calculated from the 
energy and the polar angle of the scattered electron. The pseudo-rapidity rjq of the struck quark can be 
expressed in terms of the kinematical variables x and as: 



where Ep(Ee) is the proton (electron) beam energy and s is the square of the centre of mass energy. The 
scattered electron and the struck quark are (in the QPM) back-to-back in azimuth, i.e. : 



of ±5%. 



model H. 




(5) 




(6) 



To look at deviations from these expectations, the following two variables are used: 



A?/ = 77 - Vq 

A(j) ^ (f)- (j)q 



(7) 
(8) 



where 77 (0) denotes the pseudo-rapidity (azimuthal angle) of a particle or calorimetric cluster. 
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Figure 2: Measured transverse energy flow relative to the calculated struck quark direction. Shown are non-diffractive 
("non-diflfr. DIS") and diffractive ("diflfr. DIS") data, both corrected for detector efi'ects. Also indicated are the models for 
DIS (ARIADNE) and for diffractive DIS (RAPGAP), which have been used to do the correction. Data with x < 10~^ are 
shown in (a), those with x > 10"'^ in (b). In both figures the proton direction corresponds to large positive values of Arj. 
In (b) the meaning of the 4 estimators for the energy flow is also indicated using the non-diffractive data. 

The measured transverse energy flow (integrated over azimuthal angle) around the expected 
nai've QPM direction {Ar/ = 0) is shown in figure ^ Displayed are the energy flows for non-diffractive 
as well as for diffractive DIS, in two regions of a; {x < 10~^ and x > 10"'^). The measured transverse 
energy flow relative to the nai've QPM prediction can be separated into that around Ai] — (this will 
be denoted the "current region" in the following) and the remainder of the energy flow at positive values 
of Arj towards the direction of the proton remnant (the "forward region" ) . The properties of the energy 
flow depend strongly on the kinematics, as may be seen in figure ^. The maximum of the transverse 
energy in the current region is shifted to positive values of Arj for all DIS data. The measured shape of 
the energy flows for the diffractive and the non-diffractive case is found to be very similar in the current 
region, whereas in the forward region a reduced amount of transverse energy is expected for diffractive ep 
DIS with a leading colourless remnant compared to the non-diffractive process. Also shown are the two 
models used to correct for detector effects, ARIADNE (version 4.03) for non-diffractive and RAPGAP 
(version 1.3) for diffractive DIS. 

To investigate in more detail the dependence of the hadronic flnal state on kinematical variables and to 
discuss the observed similarities as well as the differences for the two cases, 4 estimators of characteristic 
properties of the measured transverse energy flow are deflned as illustrated qualitatively in flgure ^ (b) . 
These estimators are calculated for each event using the measured calorimetric clusters and correction 
factors in case of the data and stable particles in the case of the model calculations. 



Firstly, the deviation At^^^ in pseudo-rapidity of the maximum in transverse energy from the expected 
nai've QPM direction is determined. Secondly the magnitude Expeak of the energy flow at this position is 
calculated. Next the width (Tpeak in pseudo-rapidity of the energy flow around the maximum is quantified. 
Finally the level of transverse energy Exforw away from the current region towards the direction of the 
proton remnant is determined. 

A77g^ is calculated as an i?T-weighted average of Arj in the region with \Ari\ < 2 around the expected 
quark direction, restricting the range in A^ to values with |A(/)| < 1.5: 

|Ar,|<2,|A0|<1.5 

A77e^ = ^ (9) 

|Arj|<2,lA0|<1.5 



Having calculated the observed deviation for an event, the magnitude Expeak in transverse energy 
at this position (normalized to one unit of pseudo-rapidity) is determined by adding up the transverse 
energy in a region of ±0.25 in pseudo-rapidity: 

Expeak ^ Et (10) 

|Ar;-A?)^^|<0.25 



The width apeak in pseudo-rapidity of the energy flow in the current region is obtained by determining 
the r.m.s of the Arj distribution weighted with Et around the measured position Arj^^ of the maximmn 
in a Ar] range of ±1: 



E Et ^ (A.) - A,j,„)' 

Ar;-A?7^„|<l 

E E^ 

|Ar;-Ar;^^|<l 



(11) 



Erf orw (normalized to one unit of pseudo-rapidity) is determined in the region starting one unit of 
pseudo-rapidity forward of the measured position of the maximum Et in the current region (77, + Ar^^^) 
up to a fixed value of pseudo-rapidity (77 = 3) - to stay away from the acceptance limit of the LAr 
calorimeter. 

E 

(jj,+A?)^„ + l)<»?<3 



ETforw — 7 



1 



3 - (77, + A77,„ + 1) 



Et 



(12) 



7 Results 



In this section, the dependence of the measured transverse energy flow on kinematical variables is stud- 
ied for diffractive and non-diffractive DIS, using the 4 estimators defined in the previous section. First 
the X dependence of the estimators for diffractive and non-diffractive data is compared. Next, the 
non-diffractive measurements are confronted with the most recent versions of models for DIS. The char- 
acteristics of the hadronic final state in diffractive DIS are then investigated for evidence of significant 
discrepancy from the expectation given by a partonic process with the related effects of perturbative 
QCD. To this end, first a phcnomcnological model of ep DIS in the same kinematic region is investi- 
gated, and then the diffractive final state is compared with a set of different models of deep-inelastic ep 
diffraction. 

Figure ^ shows the dependence of the estimators on x for diffractive and non-diffractive DIS. The 
measured values, together with the statistical and systematic errors, can be found in the appendix in 
tables |l| and ^. For both cases a significant deviation Arjev of the maximum in transverse energy from the 
naive QPM expectation is observed which strongly increases with decreasing values of x. For diffractive 




Figure 3: Measured estimators of the transverse energy fiow for non-diffraetive ("non-diffr. DIS") and diffractive ("diffr. 
DIS") DIS, compared with the ARIADNE and RAPGAP models. Shown is the dependence on x for the measured deviation 
from the calculated direction (a) , for the magnitude of the current region (b) , for the width of the current region (c) and 
for the forward transverse energy (d). 

DIS the measured value is found to be smaller by about 0.2 units in pseudo-rapidity than in the case of 
non-diffractive DIS. Expeak increases with increasing values of x with no difference visible between the 
two classes. No x-dependence is observed for (jpeak which has a value of w 0.4 (this corresponds to a 
full-width-half-maximum for a Gaussian shaped distribution of about 0.9). Significant differences arc 
observed in Exforw, where for non diffractive data Etjotw ~ 2 GcV/unit-of-rapidity with almost no 
dependence on x. For the diffractive case Exforw increases with increasing values of x, the magnitude 
being lower by 20 - 60 % compared with non-diffractive DIS, as expected and demonstrated clearly in 
the following. Also shown are the predictions of the models (ARIADNE 4.03 and RAPGAP 1.3) that 
have been used to correct for detector effects. The data are described reasonably well for both cases. 
At small values of x (< 10~^) the ARIADNE model (version 4.03) overestimates the measured values of 

Expeak ■ 




Figure 4: Measured dependence of the estimators on x for non-diffractive DIS, compared with the prediction of the 
ARIADNE 4.07, HERWIG 5.8d and LEPTO 6.3 models. Shown are the measured deviation from the calculated direction 
(a), the magnitude of the current region (b), the width of the current region (c) and the forward transverse energy (d). In 
addition a QPM calculation including hadronization is shown in (a). 

Decreasing values of x correspond on average to an increase in the invariant mass W of the hadronic 
system. This is expected to lead to an increase of the phase space for effects of perturbative QCD 
and particle production in the final state. The measured deviation from the naive QPM expectation is 
sensitive to this increase and also to details of the implementation of QCD effects in the different models. 
In a previous analysis it has been shown that the level of transverse energy in the central region of 
the 7*p system (this corresponds to the forward region in the laboratory frame) increases with decreasing 
X for constant Q^. 

Before investigating in more detail the diffractive contribution, some recent versions of models for 
DIS will be compared with non-diffractive data. The development of these versions of the models was 
motivated by the previously unsatisfactory description of the measured hadronic final state for x < 10~^ 
as shown in[^l|, The comparison of the predictions of ARIADNE (version 4.07), HERWIG (version 



5.8d) and LEPTO (version 6.3) with the measured data in figure ^ shows clear deviations, the differences 
being most pronounced in the estimators for Ayye^ and Exforw The shape of the dependence of Aij^v 
is very similar for these three models but not as steep at small x as in the data. The predicted shape 
and the magnitude of the Exforw dependence on x differs between models. The description of the data 
obtained with these models is in general worse than that given by version (4.03) of ARIADNE (figure ||). 
In figure ^ (a) also a QPM like calculation for the final state is shown. For this calculation (done with 
the LEPTO model) only the contribution from the Born term for DIS was considered together with 
hadronization as implemented in the Lund string model. For all values of x this calculation significantly 
underestimates the measured deviation. It should be noted that a pure Born term calculation at the 
parton level (i.e. no hadronization) gives Ayyet, = 0. This underlines the well known need to include effects 
of perturbative QCD (i.e. emission of gluons) in the modelling of the hadronic final state. However the 
understanding of the final state in the new kinematic domain of DIS opened by HERA [x < 10^^) still 
remains a challenge. 



> 1 
^ 0.9 

0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




a) 



0.2 0.4 0.6 0.8 





2 


Cl 


1.8 


D 


1.6 


i_ 


-»-' 


1.4 


'c 


1.2 


Z5 


> 


1 




0.8 


o 


0.6 


O 


0.4 


s 
o 


0.2 


1- 

II , 1 






b) 



□ HI diffr. DIS 
- MINIHERA 



□ 

' r -I- -i - 



0.2 0.4 0.6 0.8 



IS 



Figure 5: Measured estimators for diffractive DIS, compared with a model calculation for ep collisions with a reduced 
proton beam energy. The dashed lines indicate the range of the MINIHERA prediction for energies between 0.82 GeV and 
8.2 GeV. Shown is the dependence on (3 [x for the MINIHERA model) for the measured deviation from the calculated 
direction (a) and for the forward transverse energy (b). 

The properties of diffractive DIS are now investigated in more detail. As the phase space for hadron 
production depends on the invariant mass of the final state, is expected that a large part of the differences 
in the Et flow observed in the diffractive and non-diffractive case are due to the fact that in diffractive 
DIS there is an isolated leading colourless remnant. This expectation is tested by comparing a model 
calculation (called "MINIHERA") for deep-inelastic scattering of 26.7 GeV electrons and / • 820 GeV 
protons to the diffractive data. For a value of / = 0.003 the average values of W^mj^ihera ^'^d i^^ 
the diffractive data) are about equal (« 80 GeV^). This corresponds to the average value of xp in the 
data (i.e. the average momentum of the pomcron using the picture of figure |l|). To estimate the effect of 
the spread in xp as in the measurements, the calculation was done at the values / ~ 0.001 and / = 0.01. 
For all the calculations the ARIADNE model (version 4.03) was used with all other parameters identical 
to those used in the calculations for the nominal HERA conditions. The results were cross-checked by a 
calculation with LEPTO (version 6.1) which leads to the same conclusions. 

Figure ^ shows the /3 dependence of the two estimators Arj^y and Exforw for the diffractive data. 



Both estimators show a strong increase with decreasing values of f3 which corresponds to x in DIS. 
Low values of /3 correspond to large masses Mx of the ^*1P system. The limit (3^1 corresponds to 
Mx — > which means that the available phase space vanishes completely. The dependence of Arjev on 
f] in the diffractive data is described reasonably well by the MINIHERA model calculation, whereas the 
Erforw dependence is significantly underestimated for values of /? < 0.2. This calculation is used only 
to demonstrate the effect on the production of hadrons given the different invariant mass of the final 
state in the diffractive data when compared to the non-diffractive data. The obvious limitations of the 
calculation are the fact that the target is simply chosen to be a proton and furthermore a fixed value 
of / (corresponding to xp) is used, in contrast to the xp distribution in the data. However it can be 
concluded that the differences observed between the two classes of DIS data (figure ||) can be interpreted 
as a result of the differences in the available phase space for hadron production in the final state. More 
meaningful quantitative comparisons of a model of diffraction are made in the following. In this model 
the incident electron probes a pomcron with a more realistic partonic content, and therefore some of 
these shortfalls do not occur. 

In figure [g] the dependence of Ary^^ and Et form on xp is shown together with the dependence on (3 
(as shown in figure |^) for diffractive DIS. The corresponding values are listed in tables ^ and ^ in the 
appendix. The measured deviation Arj^v is observed to increase moderately whereas Exforw strongly 
increases with increasing values of xp , as is to be expected as the isolation of the colourless remnant in 
diffraction from other hadrons is reduced. 

Also shown are the predictions of three models which can generate a diffractive contribution to DIS. 
RAPGAP (version 1.3) is able to give a good description of the data over the whole range in xp and 
(3. Here diffractive DIS is modelled as deep-inelastic electron scattering of a partonic constituent of a 
pomeron, the latter being emitted from the proton. In the evolution of the final state, perturbative QCD 
effects are taken into account. A QPM like calculation is also shown. This calculation was performed 
with the RAPGAP model, considering only the Born term (i.e. electron-quark scattering) and subse- 
quent fragmentation via the Lund string model. This calculation is observed to underestimate the data 
considerably in kinematic regions (/3 < 0.3 or xp > 0.01) which corresponds to larger values of Mx and 
thus indicates the need to include effects of perturbative QCD for diffractive DIS as well. Compared 
to RAPGAP the other two models give a worse description of the data, reproducing neither the shape 
nor the magnitude of the measurements correctly. It should be noted that these discrepancies can be 
correlated with the unsatisfactory description of the non-diffractive final state in these models. 

A previous analysis |l^ by the ZEUS collaboration comparing the energy flow -j^p^ relative to the 
naive QPM expectation found at most a small deviation from this expectation for events with a LRG and 
a maximum value of the deviation of about 0.4 in Arj for events without a gap. In the analysis presented 
here signiflcant deviations from the QPM expectation are observed for both classes of DIS events. This is 
due to the choice of Et as the weight for the pseudo-rapidity distribution relative to the expected QPM 
direction. Using Et rather than E as weighting factor leads to a greater sensitivity to deviations from 
the QPM expectation for low values of x. 

As demonstrated above, the reduction of the phase space available for the final state in the case of 
diffractive DIS is the main reason for the observed differences compared with the non-diffractive case. 
This conclusion was recently also reached by the ZEUS collaboration in||l^ , where charged particle spectra 
in DIS have been analyzed in the current region of the 7*p centre of mass system. A comparison of the 
measured transverse momenta with data from fixed target experiments at a value of W comparable to 
Mx in the LRG events gave good agreement. 




Figure 6: Measured estimators for diffractivc DIS compared with several different models. Shown is the dependence on 
xp for the measured deviation from the calculated direction (a) and for the forward transverse energy (b) as well as the 
dependence on /3 for the deviation (c) and for the forward transverse energy (d). In addition a QPM calculation including 
hadronization is shown. 

8 Summary and Conclusions 

Measurements of transverse energy flow Et in deep-inelastic ep scattering have been made using data 
taken at HERA with the HI experiment. The energy flow was analyzed in the laboratory frame of 
reference for diffractive and non-diffractive data. The diffractive data are selected experimentally by 
demanding a large rapidity gap in the hadronic final state around the proton remnant direction, making 
measurements of diffraction possible in the range x < xp < 0.02. Estimators which quantify features of 
the topology of the Et flow, corrected for detector effects, have been compared with the expectations of 
different models based on QCD. 

The measurements indicate that the interpretation of deep-inelastic scattering as the scattering of a 
current quark with associated effects of perturbative QCD continues to be valid for the hadronic final 



state of the difFractivc process. The level at which these effects occur is consistent with the reduced phase 
space available in the diffractive process compared to that in non-diffractive DIS. 

The measured Et flow for diffractive DIS is well described by a model (RAPGAP), in which the proton 
couples at low momentum transfer squared t to a colourless object (pomeron). Here the deep-inelastic 
scattering process involves the partonic structure of the pomeron. 

Models for deep- inelastic ep scattering (LEPTO and HERWIG) in which the interaction of the electron 
involves the partonic structure of the proton, and not the one of an entity such as the pomeron, have 
been investigated. Here the diffractive configuration occurs because of non-perturbative QCD effects in 
the formation of the final state. These models do not describe the measurements of Et flow in diffractive 
DIS as well as RAPGAP. The observed discrepancies are however at a level which is similar to the 
disagreement observed for the bulk of DIS data, and therefore it is possible that further developments in 
these models may rectify this disagreement. 
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Appendix 



logiox 


Ariev [rapidity] 


Expeak [GeV/ rapidity] 


(Jpeak[rapidity] 


Etj 


orw [GeV/ rapidity] 


-3.72 


1.05 ±0.01 t°om 


3.5 ±0.1 1^3 


0.41 ±0.01 


1 


96 ±0.04 toio 


-3.57 


0.95 ±0.01 


3.9 ±0.1 tg:3 


0.41 ±0.01 IHi 


1 


92 ±0.03 ig:io 


-3.42 


0.84 ±0.01 ±0.02 


4.2 ±0.1 ±0.3 


0.40 ±0.01 ±0.01 


1 


94 ±0.03 ±0.10 


-3.26 


0.81 ±0.01 ±0.03 


4.6 ±0.1 ±0.2 


0.43 ±0.01 ±0.01 


1 


94 ±0.03 ±0.09 


-3.11 


0.75 ±0.01 ±0.03 


4.8 ±0.1 ±0.3 


0.40 ±0.01 ±0.01 


1 


89 ±0.03 ±0.10 


-2.96 


0.75 ±0.01 ±0.05 


5.2 ±0.1 ±0.3 


0.42 ±0.01 ±0.01 


1 


88 ±0.03 ±0.10 


-2.80 


0.71 ± 0.01 ± 0.06 


5.5 ±0.1 ±0.3 


0.42 ±0.01 ±0.01 


1 


92 ±0.04 ±0.10 


-2.65 


0.66 ±0.01 ±0.08 


5.6 ±0.2 ±0.4 


0.42 ±0.01 ±0.01 


1 


99 ±0.05 ±0.10 


-2.50 


0.63 ±0.01 ±0.08 


5.7 ±0.2 ±0.4 


0.43 ±0.01 ±0.01 


1 


83 ±0.04 ±0.10 


-2.34 


0.60 ±0.01 ±0.10 


6.6 ±0.2 ±0.6 


0.42 ±0.01 ±0.02 


1 


82 ±0.06 ±0.09 


-2.19 


0.54 ±0.02 ±0.11 


6.9 ±0.3 ±0.4 


0.41 ± 0.01 ± 0.01 


1 


74 ±0.07 ±0.10 


-2.04 


0.53 ±0.02 ±0.13 


7.6 ±0.4 ±0.6 


0.41 ± 0.01 ± 0.02 


1 


93 ±0.12 ±0.18 



Table 1: Values of the energy flow estimators as a function of x for non-diffractive DIS. The first error 
given is the statistical, the second is the systematic error. 



logwx 


A?7et, [rapidity] 


Erpeak [GeV/ rapidity] 


c^peak [rapidity] 


Et forw [GeV/ rapidity] 


-3.64 


0.77 ±0.02 ±0.05 


3.2 ± 0.2 ± 1.0 


0.38 ±0.01 ±0.01 


0.79 ±0.06 ±0.07 


-3.31 


0.62 ± 0.02 ± 0.06 


4.0 ±0.2 ±0.6 


0.41 ± 0.01 ± 0.01 


0.78 ± 0.05 ± 0.06 


-2.98 


0.56 ±0.02 ±0.07 


4.9 ±0.2 ±0.3 


0.37 ±0.01 ±0.01 


0.90 ±0.05 ±0.08 


-2.65 


0.48 ±0.02 ±0.12 


5.6 ±0.3 ±0.4 


0.37 ±0.01 ±0.01 


0.96 ±0.07 ±0.11 


-2.32 


0.46 ±0.03 ±0.18 


6.6 ±0.4 ±0.5 


0.37 ±0.01 ±0.02 


1.10 ±0.11 ±0.06 


-1.99 


0.48 ±0.09 ±0.29 


7.5 ±0.7 ±0.8 


0.35 ±0.02 ±0.03 


1.43 ±0.30 ±0.13 



Table 2: Values of the energy flow estimators as a function of x for diffractive DIS. The first error given 
is the statistical, the second is the systematic error. 





A;/, [i-aj)/(l/l jj] 


Er forw [G( 1 '/ rapidity] 


-2.86 


0.43 ±0.01 ±0.05 


0.12 ±0.01 ±0.02 


-2.57 


0.48 ±0.01 ±0.08 


0.22 ±0.01 ±0.02 


-2.29 


0.56 ±0.02 ±0.10 


0.46 ±0.02 ±0.03 


-2.00 


0.61 ±0.02 ±0.15 


0.77± 0.03 ±0.04 


-1.71 


0.57 ±0.02 ±0.08 


1.25 ±0.05 ±0.07 


-1.43 


0.70 ±0.03 ±0.12 


1.77±0.09±0.11 


-1.14 


0.80 ±0.05 ±0.15 


2.52 ±0.17 ±0.66 



Table 3: Values of the energy flow estimators as a function of xp for diffractive DIS. The first error 
given is the statistical, the second is the systematic error. 



/3 


Arjev [rapidity] 


Et forw [GeV /rapidity] 


0.07 


0.79 ±0.01 ±0.09 


1.54 ±0.04 ±0.26 


0.21 


0.58 ±0.01 ±0.07 


0.61 ±0.02 ±0.05 


0.36 


0.45 ± 0.01 ± 0.06 


0.32 ±0.01 ±0.01 


0.50 


0.34 ±0.01 ±0.06 


0.14 ±0.01 ±0.01 


0.64 


0.23 ±0.01 ±0.04 


0.09 ±0.01 ±0.01 


0.79 


0.15 ±0.01 ±0.03 


0.04 ±0.01 ±0.01 


0.93 


0.05 ±0.01 ±0.01 


0.01 ±0.01 ±0.01 



Table 4: Values of the energy flow estimators as a function of (3 for diffractive DIS. The first error given 
is the statistical, the second is the systematic error. 
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